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Abstract. Motivated by recent experiments [Claussen N R et al. 2003 Preprint 
cond-mat/0302195 we investigate the magnetic-field dependence of the Josephson 
frequency of coherent atom-molecule oscillations near a Feshbach resonance. Far off 
resonance this frequency is determined by two-body physics. However, the Josephson 
frequency is relatively large and a description purely in terms of scattering lengths turns 
out to be inadequate. In particular, we have to include the effective range parameter 
of the interatomic interaction in our calculations. Close to resonance, the frequency 
deviates from the two-body results due to various many-body shifts. Considering also 
the many-body effects, we find perfect agreement with the experimental results over 
the entire range of magnetic field. 
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1. Introduction 

The collisional properties of ultracold alkali atoms are strongly affected by the presence 
of a bound state in a coupled closed channel when the energy of that bound state is 
close to the energy of the two colliding atoms. This situation is called a Feshbach 
resonance [TJ El E] and is currently of great experimental interest, because the energy 
difference between the two colliding atoms and the molecular bound state is magnetic- 
field dependent due to the different Zeeman shifts of the atoms in the open channel and 
the molecule. This leads to a high level of experimental control over the interatomic 
interactions j3j. 

The physics of a Feshbach resonance bears some similarity with an ordinary 
quantum-mechanical two-level system, where one level represents the two atoms and the 
other level the molecular state. Inspired by this picture Donley et al. have conducted 
a Ramsey- type experiment, by starting from an essentially pure atomic condensate and 
applying two short pulses in the magnetic field separated by a longer evolution time [Sj. 
As a function of the evolution time oscillations in the number of remaining condensate 
atoms are observed. Sufficiently far from resonance the frequency corresponds to the 
molecular binding energy, in agreement with the fact that in that case the coupling may 
be neglected and the Josephson frequency of the oscillations between an atomic and a 
molecular Bose-Einstein condensate is equal to the energy difference between the two 
atoms and the molecule. More recent high-precision experiments by Claussen et al, over 
a larger range of magnetic fields, have shown that close to the resonance the frequency 
deviates from the two-body result as a result of many-body effects [6 . 

In order to be able to theoretically describe the magnetic-field dependence of 
the observed frequency of coherent atom-molecule oscillations over a large range of 
magnetic field, we have to incorporate into the many-body theory the relevant two- 
atom physics exactly, and, in particular therefore, the correct molecular binding energy. 
Far off resonance this binding energy is relatively large and it turns out to be necessary 
to include the energy dependence of the interactions and coupling constants in the 
theoretical description. In particular, the use of solely the scattering length is inadequate 
and we also need to include the effective range of the interatomic interaction potential. 
In addition, to describe the deviation from the two-body result close to resonance, 
the theory has to incorporate the correct mean-field interactions between condensate 
and noncondensate atoms that ultimately give rise to mean-field shifts of the molecule. 
Although a large number of important contributions have been made in formulating 
a theory for Feshbach-resonant interactions [3 El 13 UH EH EH EH), it appears that 
none of the existing theories provides a quantitative agreement with the experimentally 
observed Josephson frequency over the full range of the magnetic field. It is the aim of 
this paper to put forward such a theory and thereby also improve upon our previous 
calculations which were restricted to the region close to resonance |14j . 

To this end, the paper is organized as follows. In section [2] we first present the 
many-body theory that incorporates the molecular binding energy exactly. This theory 
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includes the effective range of the interatomic interactions, which is determined from 
experiment in section |3] by calculating the molecular binding energy in vacuum. In 
section E] we calculate the frequency of the coherent atom- molecule oscillations, i.e., the 
Josephson frequency that includes the many-body shifts of the two-body results. We 
end in section with our conclusions. 



2. Atom-molecule coherence 



In this section we present the theory for the description of Feshbach- resonant interactions 
between atoms. The theory presented here improves upon our results presented in ^3] 
by including the energy dependence of the various coupling constants and, in particular 
therefore, the effective range of the atom-atom interactions. 



2.1. Microscopic theory 

Our starting point is the microscopic atom-molecule hamiltonian for the description 
of a Feshbach resonance ^3]. For a homogeneous system in a box of volume V this 
hamiltonian reads 

H = Ck4«k + ^7 E V ( k - k ')aL/2+k a K/2-kaK/2+k'aK/2-k' 
k ZV K,k,k' 

+ -77 + S b{B) blb k + -JfjY.9B(k) [&K«K/2+kflK/2-k + h.c.j . (1) 



K,k 

Here, the operator a k creates an atom with momentum Kk, and its hermitian conjugate 
(h.c.) annihilates an atom with this momentum. The single-atom dispersion is given 
by ek = fr 2 k 2 /2m. The molecules are described by the operators b\ and and have 
a single-particle dispersion given by ey_/2 + 5b{B) since their mass is twice the atomic 
mass m. Here Sb{B) is the detuning of the bare molecules, i.e., the energy difference 
between the bare molecules and the atoms, which depends on the magnetic field B due 
to the fact that the magnetic moment of the molecule is different from that of the atoms. 
The Fourier transform of the atom-atom interaction in the incoming channel is given by 
V^(k), and goes to zero as k — > oo due to the nonzero range of the atomic interactions. 
The bare atom-molecule coupling is given by <?B(k), and also goes to zero as k — ► oo, 
due to the fact that the wave function of the bare molecule has a finite, nonzero extent. 
Both the atom-molecule interactions, as well as the molecule-molecule interactions have 
been neglected, because we are interested in the Josephson oscillations of an almost pure 
atomic condensate. 

In applying the above microscopic hamiltonian to realistic atomic gases we have 
to perform perturbation theory in the atom-atom interaction and the atom-molecule 
coupling. These quantities are large and we have to consider them to all orders. At 
the low densities of interest to us here, we need only to consider all ladder diagrams, 
since three and more body processes are negligible. This automatically builds into our 
theory the relevant two-atom physics, i.e., the scattering amplitude of the atoms and the 
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Figure 1. (a) Born scries for the scattering amplitude. (b) Diagrammatic 
representation of the many-body T matrix. The solid lines correspond to single-atom 
propagators. The dashed lines corresponds to the interatomic interaction ^(k). 



binding energy of the molecule. Note that the fact that both the atom-molecule coupling 
and the atom-atom interaction decrease rapidly with increasing momenta ensures that 
no ultraviolet divergencies are encountered when doing perturbation theory in these 
quantities. In the next section we discuss the required perturbative expansion. 



2.2. Ladder diagrams 

First, we consider the atom- atom scattering amplitude. To lowest order in the 
interaction the amplitude of two atoms scattering from relative momentum hk.' to 
relative momentum frk is simply equal to V(k — k'). However, since this interaction 
is strong we have to consider also the situation where the atoms collide twice, three 
times and so on. Diagrammatically, this leads to the Born series shown in figure Q (a), 
which is summed by introducing the many-body T(ransition) matrix that obeys the 
Bethe-Salpeter equation 

T MB (k, k', K, z) = V(k-k') 



1 + iV(e K / 2 +k" - A*) + iV(eK/2-k" - n) 



k" 

xT MB, 



Z — 6K/2+k" — e K/2-k" 

! (k",k',K,^) . (2) 

Here we have used the Hartree-Fock approximation for the thermal atoms and, for the 
moment, neglected mean-field shifts. Furthermore, N(x) = [e^ x — is the Bose 
distribution function of the atoms, /i their chemical potential, and f3 = l/k^T the 
inverse thermal energy. The total energy at which the scattering takes place is denoted 
by z. Due to the surrounding medium the scattering amplitude depends also on the 
center-of-mass momentum K of the colliding atoms. The diagrammatic representation 
of this equation is given in figure [T] (b). 

For temperatures not too close to the critical temperature we are allowed to neglect 
the many-body effects JTS], and (J2J) reduces to the Lippmann-Schwinger equation for 
the two-body T matrix. The effective scattering amplitude thus becomes T 2B (k, k', z — 
ck/2). For the small external momenta of interest to us here the momentum dependence 
of the two-body T matrix may be neglected. Its energy dependence, however, cannot be 
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Figure 2. Diagrammatic calculation of the amplitude of two atoms to form a 
molecule. The open and filled triangles represent the bare and many-body coupling 
constant between the atoms and the molecule, respectively. 



neglected, since we are ultimately interested in calculating the molecular binding energy 
which is relatively large far off resonance. Including the energy dependence of the T 
matrix we thus conclude that the renormalization of the amplitude for two atoms 
scattering at the physically relevant energy hu + = hu + iO is given by 



V(k-k') ^T 2B (0,0,^ H 



4:7iab K h 2 



in 



1 



1 + lOhe 



2h 



(3) 



Here a^g denotes the s-wave scattering length of the potential V(k), which, due to 
the fact that it describes the nonresonant part of the atomic interactions, is called the 
background scattering length. The parameter r^ g is, by definition, the effective range of 
this potential. 

Next, we calculate the amplitude for two atoms with momenta HK./2 + hk and 
KK./2 — hk to form a molecule with energy z. To lowest order in the atomic interaction 
this amplitude is equal to g B (k). However, before forming a molecule, the two atoms 
may undergo multiple collisions and we again have to take into account the full Born 
series, diagrammatically given in figure El (a), which is summed by again introducing 
the many-body T matrix as shown in figure El (b). This leads to the equation 

1 



9 



MB, 



k,K,z)=g B (k) + ^T 
v k' 



UB (k,k',K,z) 
1 + iV(e K / 2 +k' - A*) + iV(e K / 2 -k' 



/i) 



<?B(k') 
Z ~ CK/2+k' — CK/2-k' 

Neglecting many-body effects, the amplitude becomes g 2B (k, z — ck/2) with 



g 2B (k,z)=g B (k) + v Y,T m (k,k>,z) 



z-2ei 



-9bW 



(4) 



(5) 



From the above equation we infer that the energy dependence of the amplitude is the 
same as that of the two-body T matrix. This result is easily understood by noting that 
for a contact potential V(k) = Vq and we simply have that g 2B = g B T 2B /Vq. Hence we 



Josephson frequency of resonantly- coupled atomic and molecular condensates 



6 



have for the renormalization of the amplitude for forming a molecule with energy %u 
the expression 



1 



l + m I mm. _ 



(6) 



g B ^ g 2B (0,huj + ) =g 

I - in, 

h 2h 

where g is the atom-molecule coupling constant at zero energy. 

Finally, we have to take into account the fact that the coupling between the atoms 
and the molecule gives the molecule a finite lifetime and leads to a change in the bound- 
state energy of the molecule. Physically, this is described by the self-energy of the 
molecules which, within our approximations, is equal to TiY^(z — ck/2) with 

' /k <A0,2e k ) 2 -4r- (7) 



(zTiy z — 2ek 

This expression is easiest understood by considering the fact that evaluated at positive 
energy hu + and using Im[l/(a> ± — a/)] = ^fn5(uj — a/), it immediately leads to the 
correct Fermi's Golden Rule for the decay rate of the molecule. Using the result for 
the amplitude in © the retarded self-energy hE^>(huj), i.e., the selfenergy evaluated 
at energy hu + , becomes 



-,2 



g m 
. -I>\ 1 - 2^ 

a bg 



«abe 



(8) 



where we have for simplicity assumed that the background scattering length is negative, 
which is the case for the applications of interest here. We have also subtracted the z = 
part in the integral in (JJJ) since this shift renormalizes the bare detuning according to 
Sb(B) — > 5(B). Note that the square-root behaviour of ^SW(Sw) at small energies is 
in agreement with the expected Wigner-threshold law for the molecule to decay into the 
two-atom continuum. 

At this point we make the connection with experimentally known parameters. The 
resonance is characterized experimentally by a width AB and a position B Q . More 
precisely, the s-wave scattering length of the atoms as a function of magnetic field is 
given by 

a(B) = a bg ( 1 - -^-) . (9) 



v B-B 0/ 

Therefore, in order to reproduce the experimentally observed width of the resonance 



we have that g = h^27rab g A.BAfj,/m, since an elimination of the molecular field shows 
that a^gAB = mg 2 / (2%h 2 A/x) . We have made use of the fact that the detuning is given 
by 6(B) = Afi(B — Bo) where A/i is the difference in magnetic moment between two 
atoms and a bare molecule. 

In this paper we focus on the experiments by Claussen et al. jO]. In these 
experiments the Feshbach resonance at B = 155.041(18) G(auss) in the |/ = 2;m^ = 
—2) state of 85 Rb is used. The width of the resonance is equal to A = 10.71(2) G 
and the background scattering length is ab g = — 443ao with a the Bohr radius. The 
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Figure 3. Molecular binding energy in vacuum. The solid line shows the result of a 
calculation with rb g = 185ao, the dashed line shows the result for rb g = 0. The dotted 
line shows |e(-B)| = h 2 /ma 2 . The experimental points are taken from UJ. 



difference in magnetic moment between two atoms and the bare molecule is equal to 
A/x = — 2.23/ib with hb the Bohr magneton ^U|. In the next section we also determine 
the effective range of the atom-atom interactions by calculating the binding energy of 
the molecule and comparing the result with experiment. 

3. Molecular binding energy 

The molecular bound-state energy is determined by solving for Two in the equation 

1w = 5(B) + nY$(tiuj) . (10) 

For positive detuning this leads to an imaginary bound-state energy, in agreement with 
the fact that the molecule decays when its energy is above the threshold of the two-atom 
continuum. For negative detuning there is a bound state with energy e m {B). Close to 
resonance we are allowed to neglect the contribution from the effective range to the 
self-energy in (JSJ) and we solve (fTUj) analytically to find 

h 2 

UB) ^ — rr~RX\2 ■ ( n ) 

m[a{B)\ z 

This is indeed the correct result for the molecular bound state energy close to resonance 
M- 

Far off resonance, we are no longer allowed to neglect the effective range of the 
interatomic interaction. Because of this, the equation for the bound-state energy can 
not be solved analytically anymore, but is nevertheless easily solved numerically. In 
order to determine the effective range we compare the result with the experimental 
results of Claussen et al. [Oj. In these experiments the Josephson oscillation frequency 
for coherent atom-molecule oscillations is determined as a function of the magnetic 
field. Far off resonance this frequency is essentially equal to the molecular binding 



Josephson frequency of resonantly- coupled atomic and molecular condensates 



8 



0.25 




Figure 4. Dispersion relation of the collective modes of an atom-molecule system for 
a condensate density of n a = 2 x 10 12 cm~ 3 at a magnetic field of B = 160 G. The 
upper branch corresponds to the gapless dispersion for phonons. The lower branch 
corresponds to coherent atom-molecule oscillations. The momentum is indicated in 
units of the inverse coherence length = ^/167ran a . 



energy and therefore independent of the condensate density. By comparing our results 
with the experimental data far off resonance, we determine the effective range uniquely. 
In figure |3] we show the binding energy as a function of the magnetic field. The solid 
line is the result of a calculation with r bg = 185a . The dashed line shows the result 
for zero effective range and the dotted line shows the magnetic-field dependence given 
in Qlip. The experimental points, taken from jB], are also shown. It should be noted 
that the errorbars on the experimental data are roughly indicated by the size of the 
points. Clearly, the result with zero effective range, as well as the expression given in 
(111}) deviate significantly from the experimental data for magnetic fields larger than 158 
G. Interestingly, the result for zero effective range deviates more from the experimental 
result than the expression given by ([lip. The curve obtained with an effective range of 
7"b g = 185ao agrees very well with the experimental data in the regime where B > 157 
G and hence we use this value for the effective range of the interatomic interaction 
from now on. As expected, close to the resonance, the results of the three calculations 
become identical because the energy dependence of the effective interactions becomes 
less important in this regime. Note, however, that all three results deviate from the 
experimental points for magnetic fields smaller than 157 G. This deviation is explained 
by considering many-body effects which are the topic of the next section. 



4. Josephson oscillations 

To determine the Josephson frequency of coherent oscillations between the atomic and 
the molecular condensate we calculate the excitation spectrum of the system by means 
of linear response theory. Therefore, we write for the atomic and molecular Heisenberg 
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creation operators 
a k (t) 

kit) 



e -iKt/h ; and 
—i2/xt/h 



(12) 



respectively. The operators d k and /3 k describe a Bogoliubov quasi-atom and quasi- 
molecule, respectively, and in linear response theory we keep terms up to quadratic order 
in these operators. The number of particles in the atomic and molecular condensate is 
equal to iV a and N m , respectively. In the presence of an atomic condensate it is crucial 
to take into account its mean-field effects on the thermal atoms. This mean-field energy 
is, in the Hartree-Fock approximation, given by 

/ 



+ 



HF 



2n R 



T 2B [0,0, fi- HE 



HF 



V 



g 2B (0,0,// - hE 



HP 



fr£ HF + /X - 5(B) 



hz£ ] (/i - ft£ HF ) 



(13) 



and reduces to the usual expression &na(B)Ti 2 n & /m far off resonance. Here, n a = N&/V 
denotes the density of the atomic Bose-Einstein condensate. 

To find the excitation spectrum, we have to diagonalize the hamiltonian in terms 
of the Bogoliubov quasi-particle operators. This requires the diagonalization of a 4 x 4 
matrix equation for the coherence factors w k , t> k , it k and u k . However, the coherence 
factors for the molecular operators are straightforwardly eliminated. This yields the 
2x2 eigenvalue problem given by 



Uk + 



/i + 2T 2B (2 / u-^-e k /2)n 8 



e k - fi + 2T e f(2/i + hu - e k /2)n a ] u k + 7^(2//)^ 
where the effective atom-atom interaction is given by 



Vk = - 



(14) 



T e f(M = T 2B (0, 0, hu + - 2^S HF ) 



+ 



2B 



(0, 0, Tloj~" 



— 2ft£ HF 



(15) 



Twj+ - 5(B) - hZ^ (hu - 2ft£ HF ) ' 
The chemical potential is determined by the time-independent Gross-Pitaevskii 
equation, which, for energy- dependent interactions, reads 

H = T 2B (2/,)n a . (16) 

Note that this last equation and (|13|) uniquely determine the chemical potential and the 
Hartree-Fock mean-field energy for a given atomic condensate density. Also note that it 
is crucial to include the Hartree-Fock mean-field shift to have equilibrium solutions. The 
number of molecules in the molecular condensate can be calculated from the number 
of condensate atoms, but is eliminated here to get the effective atom-atom interaction 
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Figure 5. Josephson frequency of coherent atom- molecule oscillations for an atomic 
condensate density of n a = 2 x 10 12 cm -3 . The solid line shows the result of the 
calculation of the Josephson frequency. The dashed line shows the molecular binding 
energy in vacuum. The experimental points are taken from [B]. 



The dispersion relation of the collective excitations is found by solving the 
eigenvalue problem in ()14J) numerically. This yields two branches, one gapless branch 
corresponding to the phonon modes, and a gapped branch corresponding to the atom- 
molecule oscillations. Physically, the difference between the two branches is understood 
by realizing that for the phonon modes the phases of the atomic and the molecular 
condensate are locked to each other and oscillate in phase. Since the hamiltonian is 
invariant under transformations a — > ae ld and b — > be 2td we conclude that the phonons 
are indeed gapless, and, in fact, correspond to the Goldstone mode associated with the 
breaking of the U(l) symmetry by the condensates. For the coherent atom-molecule 
oscillations the phases of the atomic and molecular condensate oscillate out of phase 
and hence the associated dispersion is gapped. Figure |U shows the two branches of the 
dispersion for a condensate density of n a = 2 x 10 12 cm" 3 at a magnetic field of B = 160 
G. 

The zero-momentum part of the gapped branch corresponds to the Josephson 
frequency of the coherent oscillations between the atomic and the molecular condensate, 
as observed by Claussen et al. . In figure El the dotted line shows the result of the 
calculation of this frequency for a condensate density of n a = 2 x 10 12 cm -3 . This density 
corresponds to the effective homogeneous density we have to take in order to compare to 
the experimental results found by Claussen et al, which are done in an inhomogeneous 
magnetic trap jE|. Roughly speaking, this comes about because in the comparison of 
our homogeneous calculation with the inhomogeneous situation the effective Josephson 
coupling, which mostly determines the deviations from the two-body result, is reduced 
by an overlap integral of the atomic and molecular wave functions. This reduction factor 
multiplied by the experimental central density determines the effective homogeneous 
density [Ti] . Clearly, our result shows perfect agreement with the experimental data 



Josephson frequency of resonantly- coupled atomic and molecular condensates 



11 



points over the entire experimentally investigated range of the magnetic field. The 
dashed line in figure |3] corresponds to the molecular binding energy in vacuum, also 
shown in figure El As expected, the Josephson frequency becomes equal to the binding 
energy far off resonance. 



5. Conclusions 

We have improved upon our previous work i n calculating the Josephson 

frequency of coherent atom-molecule oscillations. In doing so, we had to include 
the energy- dependence of the interactions, and, in particular, the effective range of 
the interatomic interactions. We have found excellent agreement with the available 
experimental results. In future work we intend to study the density and magnetic-field 
dependence of the Josephson frequency in more detail. 
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